Piezoelectric coupling, phonons, and tunneling into a quantum Hall edge 
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We show that the piezoelectric coupling to three-dimensional phonons in GaAs renormalizes the 
current-voltage exponent for tunneling of electrons into an incompressible quantum Hall edge. The 
leading correction is always negative, in agreement with experiments on the v = 1/3 state and, 
depending on the precise value of the edge plasmon speed, can be as large as a few percent. We 
also discuss higher-order corrections, which determine the effect of the piezoelectric coupling in the 
extreme infrared limit. 
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I. INTRODUCTION 

The discrepancy between the tunneling exponent pre- 
dicted by the conformal field theory^ (CFT) of a quan- 
tum Hall (QH) edge and the one measured in cleaved- 
edgc-overgrowth (CEO) experiments 2 ' 3 , for a review see 
Ref. 0, remains a puzzle. For the principal filling frac- 
tion v = 1/3, the discrepancy is not that large but is 
still believed to require an explanation (for recent dis- 
cussions of the issue, see Refs. I5l6l) . A natural way to 
resolve the discrepancy would be to identify additional 
gapless modes, not present in the original chiral CFT. 
Such modes appear, for example, in scenarios based on 
the "edge reconstruction" J£ Another group of explana- 
tions involves the role of the Coulomb interaction in the 
presence of a "hard" edgep^ili profiles of the electron den- 
sity obtained numerically in this case are rather similar 
to those resulting from edge reconstruction. Finally, al- 
though we do not address that case in the present paper, 
we note that there seems to be a discrepancy between 
theory and experiment also for compressible QH states; 
for more detail, see Refs. I ill ll 

The effect of a one-dimensional (ID) "phononlike" 
mode on scattering between two QH edges was considered 
in Ref. Il2t and it was shown that a derivative coupling of 
the charge density to such a mode does renormalize the 
resistance. For tunneling from a bulk metal into a sin- 
gle edge, as in the CEO experiments, this coupling will 
similarly renormalize the tunneling exponent. The req- 
uisite ID phononlike mode, propagating only along the 
edge, can be a result of the edge reconstruction or other 
mechanisms enumerated in Ref. Il2l all of which have to 
do with the near-edge properties of the electron density. 

For the purpose of explaining the experimental data, 
however, one does not have to consider the electronic sys- 
tem alone. In particular, in the present paper we note 
a hitherto apparently unidentified effect of the ordinary 
three-dimensional (3D) acoustic phonons in GaAs. The 
derivative coupling of the type considered in Ref. ^3 is 
irrelevant in this case, in the sense that it does not lead 
to a logarithmic correction to the propagator of the edge 
plasmon. (The replacement of the 3D phonon with a ID 
phonon in this problem in Ref. looks to us completely 
arbitrary.) As we show here, however, there is a loga- 



rithmic effect (and consequently a renormalization of the 
tunneling exponent) due to the piezoelectric coupling. 

We limit ourselves to the principal filling fractions 
v = l/(2p+ 1) and assume that the CFT contains a sin- 
gle chiral boson — the left-moving edge plasmon mode.* 
However, the phonon-induced correlation to which we at- 
tribute our results may also play a role in a much broader 
class of problems involving ID conductors. 

The piezoelectric electron-phonon interaction is de- 
scribed by the Hamiltonian 



piczo 



hd x e L 4>6(y)5(z) 



(1) 



Here 9l is the chiral field of the edge plasmon, -kzdJSi, 
is the representation of the electron density in the edge 
CFT, <j) is the phonon field, and h is the coupling con- 
stant. The delta- functions restrict the coupling to one 
dimension, and indeed Or, is a ID field; however, <p is 
fully three-dimensional. 

The piezoelectric coupling is known to be highly 
anisotropic, but here we follow a common practice and 
use a direction average. The main result — a renormaliza- 
tion of the tunneling exponent — does not depend on this 
replacement. In fact, one could get rid of the averaging 
directly in the solution, Eq. pifl below, by bringing h 2 
under the integral over k± (the transverse wavenumber 
of the phonon field) and supplying it with an angular 
dependence. This would not affect the structure of the 
infrared-sensitive terms. 

We keep in mind that there are three polarizations of 
acoustic phonons (two transverse and one longitudinal), 
but in our calculation their contributions simply add up, 
so we will think of as representing just one of them. 
Interactions of the form JQ| have been considered in a 
variety of problems concerning GaAs structures (for re- 
cent work, see Ref. Il4|) but, to our knowledge, not in 
connection with tunneling into a QH edge. 

It is straightforward to count the powers of momenta 
and see that the Hamiltonian allows for a logarithmic 
correction to the propagator of 0l. We can anticipate the 
nature of the effect by noting that the Hamiltonian gives 
rise to a non-local time-dependent "potential" between 



2 



electrons, which to the leading order in h 2 has the form 



U(x,t) = - 



) 



(2) 



Here v p h is the phonon speed, and r the Euclidean time. 
Since this "potential" is attractive, we expect that it will 
tend to confine the electron cloud near the tunneling site, 
thus leading to a decrease in the tunneling exponent a. 

What remains, then, is to compute, using the inter- 
action Hamiltonian QJ, the coefficient of the logarithm, 
to see if the effect can be large enough to be experi- 
mentally accessible. In Sect. IIIII we will see that it 
can. The precise answer depends on the ratio c_/u p h 
of the speed of the chiral mode to that of the phonon. 
For c_ in the range 10 5 -10 6 cm/s, the leading (in h 2 ) 
correction to the tunneling exponent ranges from an 
accessible Aa = —0.09 to an apparently unobservable 
Aa = —0.003. The correction is always negative, in 
agreement with the experiments^ on the v = 1/3 state. 
Higher-order corrections and the v = 1 state are dis- 
cussed in Sect. IIVI Sect. 0is a summary of the results. 



II. THE MODEL AND THE PATH INTEGRAL 

We start with the Euclidean action that consists of 
three parts: the individual actions of the ID chiral boson 
6l and the 3D phonon <f>, and the interaction correspond- 
ing to the Hamiltonian (QJ: 



The easiest way to obtain a path integral for the chiral 
field 8l is to separate it into the even and odd (with 
respect to x) components: 6l — 9 e + 6 Q . If we consider 
9 e as the canonical coordinate, then it follows from Q 
that 



p = 2nd x 9 a 



(8) 



is the corresponding canonical momentum. The path- 
integral measure can be written as T>9 e T>p or, since the 
Jacobian of transformation from p to o does not depend 
on the field, equivalently as T>9 e V9 = T>6l, Thus, the 
vacuum-to- vacuum amplitude in the presence of a current 
J(x, t) can be written as 



V9 L e 



-S E + J drdxJ9 L 



(9) 



We assume that tunneling occurs at isolated sites — 
impurities, or tunneling centers, and in what follows con- 
sider just one such site, located at x = 0. The tunneling 
term St is then 

S t = const J dr[c(Q, t)^ (0, r) + H.c] , (10) 

where the operator c(0,r) destroys a spin-polarized 
Fermi-liquid electron at point x = 0. Treating St pertur- 
batively, we can obtain the tunneling rate via a version 
of the optical theorem, from the two-point correlator of 
ip. The requisite correlator is given by Eq. jnj with a 
special choice of the current: 



St 



So + Sci, 



drd i xH n 



(3) 



A term that describes tunneling of a single electron be- 
tween a Fermi-liquid and the QH edge will be added later. 
InEq. ©, 

S e = n J dTdx[-id T 9 L d x 9 L +c^(d x 9 L ) 2 ] , (4) 

5 = y'drd 8 x|i(M) 2 + ^h(M a } > (5) 

K = 1/Attv, and t = it is the Euclidean time. We denote 
ID integrations by dx, and 3D ones by d 3 x. Our normal- 
ization of the field 9l is such that the expansion in terms 
of the canonically normalized creation and annihilation 
operators has the form 



h = y , 1 =(b k e 



■iuj—t-\-ikx _j_ ^ ^iuj-t — ikx\ 



z.m. , 



(6) 

where L is the length of the edge, and z.m. stands for the 
zero modesi In this normalization, the fcrmion creation 
operator in the chiral CFT is 



^ =■ cxp(-6» L ) : 



(7) 



J(x,t) = -[S(x-x 1 )5(t-t 1 )-6(x~x 2 )5(t-t 2 )} , (11) 

V 

corresponding to insertion of a fermion at x = sci, r = r± 
and removal of a fermion at a; = x 2 , t — t 2 . For the 
present purposes, we only need to consider x\ = x 2 = 0, 
but the structure of the correlator is elucidated by taking 
general x\ and x 2 , so we compute it for that more general 
case. 

Note that the correlator just defined corresponds to 
the bosonic version of time ordering, 

Q{X,T) = (01^(12,7^(3*, Ti)]|0) , (12) 

where x = x 2 — x\ , r = t 2 — n , and 

T B [A(T2)B( Tl )} = e(r) J 4(T 2 )B(r 1 ) + e(-r)B(r 1 )A(r 2 ) , 

,( 13 ) 

even though the exponentials J7J are in fact fermions. 
This is inconsequential since the correlator l|12l) does not 
occur as an internal line in any Feynman diagram. 

We consider the case of zero temperature, T = 0. The 
optical theorem gives the tunneling rate per a unit inter- 
val of the biasing energy E > in terms of the analytical 
continuation of G(x, r) to real time. For tunneling into 
the edge, the rate is 



^ oc N(-E)lm[i / dte lEt g{0 7 it + 5)] 
dE J Q 



(14) 
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where N(—E) is the density of states in the Fermi-liquid. 
As indicated by the infinitesimal 6 > 0, in Eq. (|14l) we 
need the values of the correlator just below the real-i 
axis. These can be obtained by analytically continuing 
the Euclidean Q(x, t) from t > 0. 

Because the theory (01 is Gaussian, the path integral 
@ can be computed exactly: 



Z[J] = exp(i JGJ) 



(15) 



where G is the full Green function of the chiral boson 9l 
[convolution integrals are implied in the exponent of Eq. 
(|15p]. In the absence of phonons, G can be replaced by 
the free Green function 



Gq(x, t) = —iv 



dQ.dK e - inT+lk * 



2tt k x (Q — ic-k x ) 
which to logarithmic accuracy equals 



G (x,t) 



v In — , t > 



V\VL , T < 



(16) 



(17) 



L is an infrared cutoff, any dependence on which will 
disappear when we use Go ( m place of G) in Eq. ltT5|l . 

At small E, where N(—E) tends to a constant, Eq. 
(|14fl then becomes 

— oc Jm[i / dte iEt (it + S)- a ] = — i?- 1 , (18) 
dE J T(a) 

where a = Ann — 1/v. The rate at a fixed biasing voltage 
V is 



eV 



dn 



Tl = I -^dE oc v a . 
lo dE 



(19) 



Thus, in the model without phonons, the tunneling ex- 
ponent is a = \/v — the prediction of the chiral CFTi 



that determines the tunneling exponent a. Consequently, 
we are interested in logarithmic contributions to G. 

Formally, Eq. H20J1 is the exact solution to our problem. 
However, it is not easy to analyze. Alternatively, it can 
be used to obtain a perturbative expansion of G, i.e., an 
expansion in powers of h 2 : 



G(x,t) =G (x,t) + G 1 (x,t) + 



(22) 



The leading term is the free propagator (|l(j|) . and the 
next-to-leading term equals 



d 3 kdn 



-iilr-\-ik x x 



An 2 J (2tt) 4 (n-ic^k x ) 2 (vl h k 2 + fl 2 ) 



Since the coupling h 2 is relatively small, we first concen- 
trate on this term; effects of the higher-order corrections 
are discussed in the next section. 
Writing 



Gi(x,r) 



h 2 



16TT 2 K 2 V 2 h 



I(x,t) 



(24) 



we find three types of logarithmic terms in the integral 
/: / k Ii + L 2 + h; for t > 



h = 
h = 



ln[A(c_r + ix) 

C-T + ix 



2 9 



ln- 



L 



h = -- 



C—T + ix 
1 



2 g t± : v P h + sc- 



ln- 



Vph'l 



isx 



A is the ultraviolet (momentum) cutoff. The approxi- 
mate sign in the expression for I means that only terms 
that give rise to lnr are kept. 

Setting x = 0, we obtain, to logarithmic accuracy 



Gx{x,t) = 



h 2 v 2 



In' 



(25) 



III. CORRECTION TO THE TUNNELING 
EXPONENT 

We are interested in corrections to the results (|18fl and 
(fT^fl due to the piezoelectric coupling Q . The full Green 
function G of the chiral boson is given by the following 
Fourier transform 



G(x,t) 



^. e -inr+ik^ mkxh (20) 



where 



F" 1 = 2hik x (itt + c_fc x ) - h 2 k\ 



d 2 k ± 1 
(2tt) 2 Q 2 + v 2 ph k 2 



(21) 

k± = (k y ,k z ), and k 2 = k 2 + k\. Note that it is the 
coefficient of lnr in the boson Green function (at x = 0) 



The corresponding correction to the tunneling exponent 
is 



Aa 



2h 2 



v ph c-{v ph + c_) 



(26) 



where the factor of 2 is due to the presence of two trans- 
verse polarizations of the phonon (the third, longitudinal, 
polarization has larger velocity and gives a smaller con- 
tribution). Note that the correction is always negative. 

For estimates, we use the same expression for the piezo- 
electric coupling and the same values of the parameters 
as m Ref . UJ In our present notation, 



1 



eei4 



Up M ^ 47TE 



(27) 



where pu = 5.36 g/cm 3 , ei 4 = 0.16 C/m 2 , and e = 
13.2eo; in this equation only, we have restored fi. With 
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these numbers, h 2 — 5.39 x 10 8 m 3 /s 3 . Using also w p h = 
3000 m/s, we can rewrite Eq. (|26H as a function of a 
single parameter, the ratio 



We obtain 



Aa 



c-A'ph 



0.04 

r(l + r) 



(28) 



(29) 



Clearly, the smaller is r the larger is the correction. 
While for c_ = 10 4 m/s it is only 0.003 in the abso- 
lute value, for c_ = 10 3 m/s it is already 0.09. The 
latter number is roughly of the same order of magnitude 
as the discrepancy between the CFT prediction a = 3 for 
v = 1/3 and the central values a = 2.65-2.85 obtained 
experimentally 2 3 for different samples. This number is 
also larger than the experimental uncertainty quoted for 
the experiment of Ref. 0- 

The most obvious way to experimentally test the 
present theory is to reduce the plasmon speed via a ca- 
pacitive coupling to external conductors. In that case, 
Eq. (|29|l predicts a further decrease in a. This predic- 
tion, however, holds only insofar as corrections of higher 
orders in h 2 can be neglected, cf. the next section. 



IV. HIGHER-ORDER CORRECTIONS 

Expanding the exact Fourier transform (|21|l in powers 
of h 2 , we obtain higher-order corrections to the plasmon 
Green function. In this section, we consider the structure 
of the perturbation series for the case c_ <C v p h, which 
is the simplest. In this case, the main contribution to 
G in a given order can be found by dropping Q, 2 in the 
integral over k±. The nth-order correction to the Green 
function at r > becomes 



In' 



A 

k x 

At x = all these terms are positive, so their sum re- 
mains the main contribution to the sum of the entire 
perturbation series. 

From now on we concentrate on x = 0. It turns out 
that the sum of (|30|l over n saturates at values of n that 
are much smaller than hq = Ac_t. For such n, the loga- 
rithm in (|30|l can be considered a slowly varying function 
of k x , so that (for n^0) 



G„(x,t) « - 
n 



v l* c - 



ln n^0 = ^n ln n^0 ^ 
n n n 



The last equality defines the dimensionless coupling £. 

It is somewhat more convenient to consider, instead of 
the sum of G„, the sum of their derivatives with respect 



to tiq. Thus, if we denote by AG the full correction to 
the Green function, AG = G — Go, then 



SAG 



[no] 

= -£eih 

71=1 



.,-1 

n 



(32) 



If n is regarded as a continuous variable, the expression 
under the sum has a maximum at n = n„, where 



n, =n e- i T i [l + O( — ) + 0(0] • 
n* 

We define a characteristic time tq by 



r 1 = c_Acxp(- 



h 2 v 



1) 



(33) 



(34) 



Eq. (|33|l can now be written simply as n» « t/tq. We 
see that tq is the infrared scale at which the perturbation 
theory breaks down. At r -c To, the sum over n is es- 
sentially discrete and is well approximated by the lowest- 
order term. On the other hand, if r is significantly larger 
than To, many terms contribute. In that case, n is indeed 
quasi-continuous. 

At even larger times, r 3> To/£, the steepest descent 
condition becomes applicable near n — n*, and using 
steepest descent we obtain 

= ^L(2nn^) 1/2 exp^r/ro + 0(n^ 2 )} ■ (35) 

The exponential growth of Eq. I35|) at large r indicates 
the presence of a plasmon state with energy 



h 2 v 



(36) 



The same conclusion can be reached by looking directly 
at the resummed expression (|21|) or, more precisely, its 
real-time version obtained by replacing Q with — ioj. (In 
the presence of a state with a negative energy, the ro- 
tation to Euclidean frequencies needs to be redefined.) 
Neglecting again the frequency dependence of the inte- 
gral over k±, we see that the interaction shifts the pole 
of F from uj = —C-k x to uj — f(k x ), where 



/(y = - c _fe(i-^A) 



(37) 



This function has extrema at k x = ±(c_ro) _1 , where it 
equals ±|-E P |. 

We interpret the presence of plasmon states with neg- 
ative energies as a reflection of the polaron effect — 
formation of a bound state of an electron and the phonon 
field. This interpretation is supported by the following 
estimate. Suppose we qualitatively describe the cumula- 
tive effect of the attractive interaction J5J over time by 
the time-independent potential 



U(x) 



U(x,T)dT — 2 



V ph\ X \ 



(38) 
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For an electron cloud with density distribution n{x) and 
spatial size of order of the "correlation length" c_r , the 
binding energy can then be estimated as 

If - /i 2 ln(Ac_T ) 

E e = - / dxdyn{x)U{x - y)n(y) ~ g • 

z J v P h c - T o 

(39) 

An electron cloud can be thought of as containing of 
order ln(Ac_To) plasmons (in the sense that this 

is how many plasmons are typically produced when an 
electron-hole pair annihilates). Dividing Eq. (|39() by this 
number, we obtain an estimate of energy per plasmon in 
agreement with Eq. l(3l)l) . 

A corollary to this argument is that any infrared effects 
associated with production of plasmons should saturate 
at time scales of order r (equivalently, length scales of 
order c_to, the size of the polaron). We therefore expect 
that in the extreme infrared, at biasing energies 

E~\E e \~ — , (40) 

the system, at any v, will cross over to the normal Fermi- 
liquid behavior. 

We now estimate the timescale tq for v — 1/3 and 
v = 1. The maximal phonon momentum A is A = 
t^D/vph, where ujd is the Debye frequency; for GaAs, 
u)d = 345 Kpi^ Given the uncertainty introduced by the 
directional averaging of h and by the imprecise knowledge 
of r = c_/t; p h, we prefer to extract the entire combination 
y = h 2 /Vp h c-, occurring in Eq. (|34|l . directly from the 
experimentally measured values of the tunneling expo- 
nent. In the limit c_ <C w P h, the lowest order correction 
(|26|l is simply Aa = — 2y. Assuming (and confirming a 
posteriori) that the perturbation theory works well for 



the v = 1/3 state, we find that, e.g., the sample^ with 
a = 2.7 corresponds to y = 0.15. For v = 1/3, Eq. (jSH) 
then gives t^ 1 ~ 10~ 9 rw£>, which is unobservably small. 
We conclude that for v = 1/3, in the experimentally rele- 
vant range of voltages and temperatures, the lowest-order 
result 1(26(1 is reliable. 

The situation is quite different for v = 1. Using the 
same value of y, we now obtain t^ 1 = 14r /ieV, close 
to the observable range of energies. The lowest-order 
result (26(1 predicts a negative correction to the tunneling 
exponent at v = 1. Experimentally, no such negative 
correction has been observed We see that this may 
be related to the breakdown of the perturbation theory 
for v = 1. 



V. CONCLUSION 

We have computed the leading correction [Eq. 1)26(1] . 
due to the piezoelectric coupling to 3D phonons, to 
the current-voltage exponent a for tunneling between 
a Fermi-liquid and a v = l/(2p + 1) QH edge. The 
correction is always negative, in agreement with the 
experiments 2 ' 3 on the v = 1/3 state, and its magnitude 
depends on the value of the edge plasmon speed. 

We have also shown that, in the experimentally rele- 
vant range of energies, higher-order corrections for v = 
1/3 are small and do not invalidate the leading-order re- 
sult, but for v = 1 they very well might. In neither 
case, however, the leading-order result (|26() represents 
the true infrared behavior: at a sufficiently low energy, 
higher-order effects will always become important. 

The author thanks A. Boyarsky, A. Chang, and V. 
Cheianov for discussions. 
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